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We present a microscopic theory of the collective modes of a "smooth" edge of a quantum Hall 
system, showing under what conditions these modes can be described as a set of independent bosons. 
We then calculate the electronic spectral function in an independent-boson model - a procedure that 
reduces to standard bosonization in the limit of "sharp" edge. The I-V tunneling characteristics 
deduced from this model exhibit, for low voltage, a power law behavior, with exponents that differ 
significantly from those of the sharp edge model. 



Understanding the character of edge excitations is cru- 
cial to the theory of the transport properties of two di- 
mensional systems in the presence of a strong perpendic- 
ular magnetic field, such as quantum Hall bars, quantum 
wires and dotsEl. An effective theory of edge excitations 
was first derived by WenH. He showed that a "sharp" 
edge (see below) is a realization of the one-dimensional 
chiral Luttingcr liquid (CLL) model, where the electronic 
spectral function exhibits, in general, a non trivial be- 
havior, leading to a density of states that vanishes, at 
low energy, as a power law. This theory has been well 
confirmed by detaiJfjd microscopic calculations^ and by 
recent experimentsa'El. The effect of the long range of 
the Coulomb interaction, which was initially ignored, has 
been recently included by Ziilicke and MacDonaldQ. 

All the above papers assumed the validity of the so- 
called "sharp edge" model, in which the density of the 
system drops sharply from the bulk value po to zero 
within a few magnetic lengths I = (Hc/eB) 1 / 2 . There 
are numerous indications that this is not always the 
correct model fof, the edge. On one hand, Hartree- 
Fock calculationsQ for strongly confined systems, pre- 
dict that, at sufficiently strong magnetic field, the edge 
undergoes a reconstruction, taking up a more extended 
shape. On the other hand, in the case of smooth con- 
finement, such as can be realized by gate electrodes, the 
electronic density is expected to have a smoothly varying 
profile (on the scale^f /), determined by classical elec- 
trostatic equilibriumtrfj. Detailed calculations using den- 
sity funclippLal theory, Thomas-Fermi theory, and other 
methodsEErtil have confirmed the theoretical validity p£ 
this "smooth edge" picture. Edge imaging experimentst£l 
have confirmed the relevancy of this description for gate- 
confined Hall bars. 

In this Communication we want to investigate the 
spectral properties of the "smooth edge" model described 
above. The mapping to a one dimensional chiral elec- 
tron liquid is not justified in this case, it, fact, a re- 
cent study by Aleiner and Glazman (AG)li3, based on 
the classical hydrodynamics approach, has shown that 



a smooth edge, in contrast to a sharp edge, supports 
multiple branches of edge waves. Of these, one is the 
usual edge magnetoplasmonllJ, and the others (infinitely 
many in the classical approach of AG) are phonon-likc 
and lower in energy than the magnetoplasmon. We shall 
show that, under the assumption of smooth density vari- 
ation, only a finite number of these phonon-like modes 
are correctly described as independent bosons. 

In order to calculate the electronic spectral function 
we rely on the strong analogy between this problem and 
that of a uniform electron gas in the partially filled low- 
est Landau level (LLL). The similarity between these two 
problems arises from the fact that in both cases, in a 
mean field approximation, the self-consistent potential is 
uniform, so that the electrons are distributed among a 
large number of degenerate orbitals at the Fermi energy. 
In a smooth edge this occurs because the nonuniform 
electronic density perfectly screens the field due to the 
external confinement potentiaMa. Recently, Johansson 
and KinaretJla have shown that a_qualitatively correct 
description of the spectral functionllj of the uniform elec- 
tron gas in the LLL at general filling factor is given by an 
independent boson model (IBM)lj in which a single lo- 
calized electron interacts with the density fluctuations of 
the system. An essentially equivalent procedure has been 
applied by Aleiner, Baranger and Glazmancil to study the 
spectral function of the two-dimensional electron liquid 
in a weak magnetic field. Finally, a formal justification 
of the IBM from diagrammatic, many-body theory has 
been provided by HaussmannEj. Encouraged by these 
successes, here we apply the independent boson model 
to the problem of the smooth edge. The resulting the- 
ory reduces to standard bosonization in the limit of a 
sharp edge, i.e., when there is only one branch of edge 
waves. When multiple branches are present, our results 
for the low energy behavior of the tunneling density of 
states are significantly different from those of the sharp 
edge model. The actual number of modes that must be 
included depends on the width of the edge, as explained 
below. 



1 



Let us begin by writing down the microscopic Hamil- 
tonian, within the lowest Landau level, in terms of den- 
sity fluctuations relative to the equilibrium density profile 
Poly)- 



-Sp(r)Sp(r')d 2 rd 2 r' 



(1) 



edge 



where the density operator (projected in the LLL) has 
been written as 



p(r) = Po{y)+Sp(r). 



(2) 



The integral in eq. (|I|) extends over the edge region which 
we take to be < x < L, < y < d, with L ^> d^> I, and 
translationally invariant along x (po(y) = for y < 0). 
The projected density fluctuation operator is given by 



8p{r) 



1 



^/nlL 



(.y-vkr + iv-vh) 
572 



(3) 



h^k 



where y^ = hi 2 , h and k are integral multiples of 2tt/L, 
ct is the creation operator of a Landau gauge orbital cen- 
tered a yu with wave vector k in the x direction. Note the 
restriction h ^ k which excludes the equilibrium com- 
ponent of the density. The kinetic energy is absent in 
eq. (0) due to projection on the LLL, and we have as- 
sumed that density fluctuations vanish in the bulk of the 
system, i.e., the bulk is incompressible. The terms lin- 
ear in Sp have vanished because of the equilibrium con- 
dition J p$(r)v{r — r')d 2 r' + V ex t{r) = constant where 
V ex t(r) is the confinement potential. Therefore the prob- 
lem is formally similar to that of a translationally invari- 
ant electron gas: the nonuniformity enters only through 
the restricted region of integration in eq. ([!]). 

The normal mode operators 5p n k are now introduced 
according to the definition 

f L dx f d 
Spnk= —e- lkx dyf nk (y)Sp(x,y), (4) 
Jo L Jo 

where f n k(y) are the solutions of the equation 



K (k\y-y>\)f nk (y>)^dyi = -^-f nk (y), (5) 
P Any- 



where Kq (y) is the modified Besscl function. They satisfy 
the orthonormality condition /„fc(y)/ m fe(2/) £fl |^dy = 
S nm , and vanish outside the interval [0, d]. Equation (|5|) 
is the eigenvalue problem solved by AG. The n-th eigen- 
function has n nodes in the y direction. In terms of the 
normal modes, the hamiltonian (Q) takes the form 



H = 



nk>0 



nk nki 



(6) 



where the operators b nk are defined via Sp nk = 
y/ kl 2 p / Lb^, and the eigenfrequencies ui nk are given by 



kve 2 1 \ nk ir, and v — 2irl 2 p is the usual filling factor in 
the bulk. 

It remains to be determined under what conditions the 
operators b n k are good boson operators. To this end we 
substitute eq. (|^) in eq. (JiJ), noting that when n <C d/l 
the gaussian factors in the integral can be replaced by 
5- functions on the scale of variation of f n k ■ We obtain 

8 Pnk ~ e - fc2;2 /4 1 J- cl_ k/2 c h+k/2 f nk (hl 2 ). (7) 

h 

(n < d/l). 

The commutator of two density fluctuations can now 
be easily calculated to be 



[Sp n k, Sp m -k] = e k 1 /2 -^ 2 -'^2(n h - k /2 - n h+k 



*fnk(hl 2 )f mk {hl 2 ) 



12) X 



kl 2 p 



in agreement with the commutation rules for bosons. In 
arriving at equation (||) we have assumed fci <C 1, and 
we have replaced the occupation number operators by 
their ground-state expectation values n k , which amounts 
to a linearization of the equations of motion around the 
equilibrium state. For the commutator [Sp n k, Sp m -k'\, 
with k 7^ k' , we find, at the same level of approximation, 
zero. This is because the commutator in question con- 
tains terms of the form c\c h , with h ^ h! , which vanish 
upon averaging in a translationally invariant (along x) 
state. 

Having thus completed the bosonization of the hamil- 
tonian, we proceed to the calculation of th fi . sp ectral 
function within the independent boson modeHEj. The 
model describes a single electron, localized at point r, 
electrostatically coupled to density fluctuations: 



Hibm= hu; nk bl k b nk + 

nk>0 

+^(r)^{r) £ M nk (y)[bl k e^ + b r 



k e- ikx ], (9) 



nfc>0 



where the matrix element M nk (y) is given by 



w / s 2e 2 „ , N kl 2 p 
M nk (y) = i—fnk{y)\—f-, 



(10) 



and ij}'(r) is the field operator that creates an electron 
in the LLL coherent state (gaussian) orbital centered 
at r. Tiie hamiltonian (^) can be solved by standard 
methodsEj, within the one-electron Hilbcrt space. The 
fermionic Green's function is obtained as 



G>(y;0 = -*(V^W,0)) 



= (1 - u (y))exp Y 



\nk>0 



1] 



(11) 
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where vo(y) = 2irl 2 po(y), and the sum over n and k in 
the exponent is restricted by the conditions n <C d/l 
and k <C I /I, which define the regime of validity of the 
hydrodynamic approximation. This result can also be 
obtained from—direct bosonization of the electron field 
operator, as ir£J. The Fourier transform of G(y,t)/2ir is 
the spectral function A > (y,tu) and gives the local density 
of states, which controls the tunneling current from a 
point contact located at position, y into the edge. From 
eq. (|ll| ) it can be easily showno that A > (y, uj) satisfies 
the integral equation 



ujA > (y,uj)= I g(Sl)A > (y,u - Sl)dSl, (12) 
Jo 



where 



g(y,si)=J2 



M nk (yf 



■5 (SI - COnk)- 



(13) 



Eq. (12), together with the conditions A > (y,u>) = for 
lu < and J Q A > (y,uj)duj = \ — vo(y), completely deter- 
mines the spectral function. This equation further im- 
plies that, at sufficiently small w, A > (y,uj) will have a 
power-law behavior 



A(y,u)~u 9 ^- X 



(14) 



if and only if the function g(y, SI) has a finite limit for 
SI — > 0. The tunneling current J, as a function of voltage 
V, will then behave as V a ^ v, °^ for sufficiently low voltage. 
Notice that this conclusion is completely general, and 
does not depend on the specific (hydrodynamic) model 
that led to the definition of g(y,Sl) in eq. ([ll]). In the 
general case, g(y, Si) could be computed from the micro- 
scopic density-density response function x(r, f , SI) of the 
edge as follows: 



d 2 r'd 2 r"v(r~ ?)v(r - ?')Imx(?, r", SI) /SI, 

(15) 



where v(r 

,2 



has the Fourier transform v(k) = 
(2iTe 2 /k) exp (— (fc^) 2 /4)c3. An important advantage of 
this microscopic formulation is that the finite lifetime of 
the collective modes (which is assumed to be infinite in 
the hydrodynamic model) would be taken into account 
through the width of the peaks in Imx- 

The calculation of the exponent g(y, 0) is easily per- 
formed within the hydrodynamic model. ._. Neglecting 
the weak nonlinearity of the n = modecj we obtain 
ffO/.O) = En/ 3 "(y)' where 



1 



My) = -fUv) 



(16) 



Although the cutoff at n = d/l introduces an uncer- 
tainty in the evaluation of the exponent at any given 
d, we emphasize that there would be no uncertainty if 



1.5 



0.5 
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FIG. 1. Electronic spectral function A>(u}) as a function 
of lo/ujq, where u>o = i>e 2 /nl, for edges of a v = 1 QH system 
with 1, 2, 3, 4, 6 and 10 modes. Logarithmic corrections 
to the edge magnet oplasmon dispersion are neglected. The 
dependence on y and po(y) have been eliminated neglecting 
the constant 1 — vo(y) and using f3 n — \/v, u) n k = ve 2 k/nn. 
The integral over k has been cut off as explained in the text. 



one used the microscopic formula ( |l5| ) for g(y, SY). Our 
approximate hydrodynamic evaluation of the exponent 
should be in good qualitative agreement with the results 
of the more accurate microscopic calculation. 

We observe that independently of the shape of the den- 
sity profile fa (y) — 1/ v, with negligible corrections aris- 
ing from the weak nonlinearity of the dispersion of the 
n = mode. Therefore in the sharp edge limit, when 
only one branch of edge waves exists, we recover the fa- 
miliar result of Wen's theory A > (ui) ~ cj 1 ^™ 1 . For n > 1 
Pn(y) fluctuates around an average value of 1/D in a way 
dependent on the form of the equilibrium density profile 
of the edge. This can be confirmed by explicit calcula- 
tion in the special case po(y) = (2/n)p atriy/y/d consid- 
ered by AG for a gate-confined electron gas, leading to 

the result (3 n (y) = \T 2 n (yfd/(y + d)) (2 - S n0 ), where 

T n (y) is the n-th Chebyschev polynomial. We conclude 
that the exponent in eq. (|lj) increases linearly with d 
and therefore that in the limit d — > oo (limit of infinitely 
smooth edge) the tunneling density of states vanishes at 
low energy faster than any power law, that is, a "hard" 
gap develops. However, it is easy to see that the power 
law behavior of eq. ( |l4| ) only holds for to -C De 2 /dn - an 
interval that shrinks to zero for d — > oo. 

In Figure 1 we present our numerical results for the 
full electronic spectral function, calculated from eq. jl^ ) 
within the hydrodynamic model for different edge widths 
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d. In contrast to the analysis of the low-frequency behav- 
ior, this calculation depends on the detailed form of the 
eigenfunctions f n k(y) and eigenfrequencies w n k- From a 
detailed study of the solutions of the eigenvalue equation 
(|J) we have found that the / n fc's can be treated as be- 
ing independent of k and the Lu n k 's to be linear functions 
of k up to a maximum wave vector k c = n/d for which 
the wavelength along the edge equals the wavelength per- 
pendicular to the edge. For k > k c the mode dispersion 
becomes approximately constant, and the wavefunction 
f n k becomes localized near the boundaries of the edge re- 
gion, giving negligible contribution to the spectral func- 
tion. The results presented in Fig. 1 have been obtained 
using the double cutoff n < d/l and k < k c : the results 
are found to be largely independent of the details of the 
cutoff procedure. 

Figure 1 shows clearly how the low energy pseudogap 
becomes more and more pronounced with increasing d 
(and, therefore, increasing number of branches of edge 
waves). For very large d the spectral function is found to 
converge to a S function centered at u>q = ve 2 /nl, which 
coincides with the simplest estimate of the potential en- 
ergy cost for the insertion of an electron into a frozen 
liquidB. 

Our results for d — > oo are in qualitative agreement 
with those obtained in refs. |is| and [2^ for the spectral 
function of the homogeneous electron gas, except that 
the latter is found to have a finite width. This happens 
because our hydrodynamic approach is unable to give the 
gapful collective modes of the homogeneous fluid phaseo, 
and hence our spectral function does not reduce to that 
of the homogeneous phase. 

In conclusion, we have performed an independent bo- 
son model calculation of the tunneling density of states 
for a smooth edge, and we have found that it vanishes 
at low frequency as a power, with an exponent that dif- 
fers significantly from the one found pin the sharp edge 
case. Recent experiments by ChangO have apparently 
confirmed the predictions of the CLL for the exponent of 
the tunneling density of states in a sharp edge. It should 
be interesting to extend these studies to see if and how 
the exponents change as the smoothness of the edge is 
varied. 
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